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The Image Enhancement Viewer has been accomplished. The final 
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SYMBOLS 

F(x) Intensity distribution in object space 

x Object coordinate 

F(w) Amplitude distribution in frequency space 
k Linear frequency 

w Radian frequency (27rk) coordinate 

G(o>) Filter characteristic in frequency space 

cdj Filter cutoff frequency 

H( z ) Amplitude distribution in image space 
z Image coordinate 

T 2 Maximum transmission level 

T t Minimum transmission level 

c Variance of the Gaussian pulse 
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1. INTRODUCTION 


25X1 A 


The work done under the present phase of the contract consisted of two main parts: first, 
the modification of the image enhancement viewer and the testing and calibration of this instru- 
ment; second , the continuance of research, in certain areas, which had been carried out by 
while working under previous contracts for the viewer. 

A description of the modified viewer, along with the methods of testing and calibration, is 
contained in the operating manual. The viewer was put in working order to the specifications 
listed in the manual. 


The following five research areas were investigated: 


1. A nonsharp, Gaussian pulse was theoretically passed through an optical, coherent system. 
The low frequencies of the pulse were filtered with a high pass occluding filter, and the final im- 
ages were plotted. 

2. A filmed nonsharp edge was experimentally filtered using the viewer, and the image of the 
edge was traced with a scanning photometer. The results of this trace were curve fitted, nor- 
malized, and plotted. 

3. Techniques (in this study, mainly exposure addition) for increasing edge contrast and 
sharpness using the viewer were experimentally examined. 

4. The application of various filtering techniques using actual aerial photographs was studied. 

5. Finally, the use of occluding spatial filtering as a tool in the area of photographic mensu- 
ration was investigated. 


From the study of these five areas, a comprehensive view of the capabilities and limitations 
of the image enhancement viewer and occluding spatial filtering may be obtained. 
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2. THEORETICAL STUDY OF THE EFFECT OF OCCLUDING FILTERS ON 

AN UNSHARP (GAUSSIAN) PULSE 


A full description of the image enhancement viewer (Fig. 1) and its basic properties may be 
found in both the operation manual and Itek Technical Report OD-61-6. In operation, the viewer 
employs coherent illumination, permitting an analysis of the system to treat it as being linear in 
amplitude. 

Since all object films used in the system are immersed in an index matching fluid (see manual) 
there should be no shift in phase, and the object may be considered solely in terms of amplitude. 

In order to obtain a general idea of the effect of occluding filtering on amplitude objects, a 
generalized case (considered Gaussian for mathematical simplicity) was considered theoretically. 
This Gaussian pulse is tranformed into the frequency plane and there filtered by a circular, cen- 
tered occluding filter of generalized dimension. The resulting distribution is then again trans- 
formed into a filtered image, which has been computed and plotted for the range at representative 
values. 

Considering a nonsharp pulse of Gaussian form (Fig. 2) 


fi(x) = T 2 - (T 2 - T t ) exp (~cx i 2 ) 


where f t (x) is the intensity distribution of a dark-centered pulse in the object plane. Its amplitude 
spectrum in the frequency plane is given by 


F(w) = 



1 

[fi(x) ] 2 exp ( -iwx) dx 


x 

Then, [fj(x)J 2 may be rewritten in the following form 


[f t (x)l* =[Tj - 


i 

(T 2 - TO exp ( - cx 2 ) J 2 



exp (-cx 2 ) 


1 

2 


Expanding, we get 
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1 - 1 - 


exp (-cx 2 ) 


1 - 


1 <l ~ ^ ex P (" cx2 ) “ T7~k ( l -?■) exp (-2cx 2 ) 


2x4 


-- 1 j 3 - (l - —) exp (-3x 2 ) 
2x4x6 1 T 2 / 


1 x 3 x 5 
2 x 4 x 6 x 8 

1 x 3 x 5 x 7 
2 x 4x6x8x10 


(*■ - exp (-4cx 2 ) 


1 — TfT " ) exp (-5cx 2 ) - ... 
l 2 , 


Letting the binomial coefficients be 


A(1 ) = 2 


A(2) = 
A(3) = 
A(4) = 
A(5) = 


2x4 

1 x 3 
2x4x6 

1x3x5 

2 x 4 x 6 x 8 

1 x 3 x 5 x 7 
2x4x6x8x10 


Then 




1 - (1 - ~) exp ( — cx 2 ) 
1 2 


i-S 

n=l 


A(n) (1 - ~) exp (-ncx 2 ) 


Putting this in the integral 


/»oo 

F(w) = I 

1 - A(n) (l - exp (-ncx 2 ) 

00 

n=l \ i2 / 


exp (-iwx) dx 


Using the cosine transform, because of symmetry, and the aperture limits 

,x f 


F(w) = 2 VT 2 " cos wxdx - 2VT7 A(n) (l - — 

J o £1 V Tl 


w s. 


exp (-ncx 2 ) cos u>x dx 


Integrating 


F(w) = 27T 2 Xf sine (wxf) - Vt 7 ^^A(n) /l - ~\ 


VT 
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This pulse may also have a light center in the form 
f*(x) = T t + (Tj - Tj) exp (-cx 2 ) 

The expansion and analysis using this pulse is very similar to that of fj(x). 

We now obtain the image of f t (x), after passing it through a symmetric, sharp cutoff occluding 
filter (see Fig. 3) centered in the frequency plane whose characteristic is 

G(u>) = 1 (ct)j ««) 

G(o>) = 0 (0 < u> < Wj) 

The distribution in the image plane is given by the following transform 
h t (z) = I F(w) G(w) exp (icox) do) 


Again, for reasons of symmetry, we may use the cosine transform 


* / co < 


F(u) cos (wz) dw 


2VT\ 


/• OO 


X{ I sine (wxf) cos (wz) dco 


v To 


| A(n) A. - [exp (-w 2 /4nc)] cos (u>z) dw 

4i) t n=l ' i ' 


r 


2Vt; 


r«i 

I Sil 


VT l 
VtT 


sine (wx f ) cos (toz) do> - -yr— x f sine (wx f ) cos (toz) dco 

X A (n)(l 


n=l 


/nc 


J" 


It 

T 2 


exp (-w 2 /4nc) cos (wz) du> 


r w l 

*^0 


exp (— w 2 /4nc) cos (wz) dca 


Applying the trigonometric identity 

sin (coxf) cos (wz) = 1/2 sin (x { + z)oj + sin (x f - z) 

and integrating, we get 


Cl ) 
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h t (z) = 


4% 


77 77 

2 + 2 


Vt~ T 

- — - Si (Xf + z )w 1 - — Si (xf - z)^ 

77 77 


VT^ 

VF 


y A(n)^l - ^7= vFrc exp (-cnz 2 ) 




Wl 


exp (-w 2 /4nc) cos (coz) dco 


Finally, the pulse image is given by 

, Vt7 

hj(z) = VT 2 ~ [Si ( x f+ z ) w i + Si ( x f “ z ) w i] 




exp (-cnz 2 ) 


^§ A(n) ( 1 ‘5) M 




exp (-co 2 /4nc) cos (wz) do> 


Rewriting the image term in the following form 


where a = 




V2c 
m = zVllc 


0> 

— = t 

COj 


dw = Wj dt 

h t (a, m) = VT^ - [Si (x^/2c + m)a + Si (x^Vlc - m)a] 


VTF 


- VT^ A(n)^l - ^ exp (-nm 2 /2) 


VTj 

VF 


exp (-a 2 t 2 /2n) cos (amt) dt 


Since c, the pulse variance is large, then XfV5c is large and the value of Si (x£i/2c ± m) 
approaches ir/2. Under this condition, the image reduces to 
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h t (a,m) = 




exp <-nm 2 /2) 






exp (-a 2 t 2 /2n) cos (amt) dt 


At this point, it is necessary to compute the values of the integral 


a 

V27rn 



exp (-a 2 t 2 /2n) cos (amt) dt 


Upon receiving these values, it is then possible to normalize the total energy of the filtered 
cases to the unfiltered case. 

The cases computed were those where a had the values of 1/ 4, 1/2, 1, and 2, and n varied 
from 1 to 5; m was taken in steps of 0.05 from 0 to 3.0; was taken in steps of 0.01 from 0 to 
0.99; and T 2 was set at 1.0. 

The normalization procedure is the following 



I(a, m) 


which is normalized to the total energy of the unfiltered pulse A p = 1(0, m) where a = 0. Therefore, 
the normalization factor is 

A P 

N = — ^ — 

I(a,m) 

The area of the filtered image, I(a,m), is approximated by Simpson’s rule for normalization 
purposes. 

Representative results of the computations were plotted (see Fig. 4, a through h). 

The parameters relating k, the linear frequency, and c, the pulse variance for a’s varying 
from 0 to 2, are plotted in Fig. 5. This relationship is defined by 

27T 2 k 2 


For the entire range of edge contrast (variation of T t ), the smaller filter sizes give the 
expected results, i.e., increased normalized contrast for increased filter size. The larger filter 
sizes give rather unexpected and varying results for changing contrast. A full explanation for 
this phenomenon cannot be made without further study. 
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OSRAM HBO-200 
mercury arc. 

546 -millimicron wave- 
length. imaged on pinhole 
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Note: Diagram not to scale, 
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Fig. 1 — Linear, coherent optical system showing focal planes of interest 
and defining the notation and location of the Fourier transform planes 
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(b) Light-centered pulse 
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Fig. 4 — Images of nonsharp Gaussian pulses of varying minimum 
transmissive level, filtered with occluding spatial filters in a linear, 
coherent optical system 
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Normalized Image Coordinates 


Normalized Image Coordinates 


(c) T, = 0.25 


(d) T 4 = 0.30 


Fig. 4 — (Cont.) 
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= (pulse variance) x (spatial filter frequency cutoff) 
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3. EXPERIMENTAL FILTERING OF A NONSHARP EDGE 


In order to substantiate the theoretical work, an experiment was designed with conditions 
approximating the analysis. 

Instead of a pulse, a filmed nonsharp edge was used as an object. This was placed in the 
object plane of the image enhancement viewer. Though the shape of this pulse had first been 
assumed to be Gaussian, subsequent traces revealed a number of un-Gaussian characteristics 
which could well affect the resulting filtered traces. 

Because an edge and not a pulse was used, the filters placed in the object plane did not have 
to be circular (the frequency distribution being one dimensional). The occluding filters used were 
pieces of drill rod whose width was accurately known. These filters were then carefully centered 
in the frequency (or filter) plane. 

Using a scanning photometer whose slit width was minimized to the point where, at maximum 
photometer sensitivity, using the largest filter, 100 percent deflection on the recorder (5-millivolt 
range) could be reached. The edge was then traced with the photometer placed in the image plane 
with its slit parallel with the edge. First, the image of the edge with no occluding filter and the 
images of the edge with filters corresponding to a = 1.5 and a = 2.5 were traced. 

The actual traces for the three cases are shown in Fig. 6. (The two filtered cases are ex- 
panded by a factor of 10 in the intensity scale.) 

The resulting filtered cases were normalized to the area under the unfiltered curve in the 
same manner as used in the theoretical traces. The resulting curves are plotted in Fig. 7. Be- 
cause the normalized cases are meant to be compared with the theoretical Gaussian, only the part 
of the edge corresponding to a pulse (center 1/V2c distance from the inflection point) is taken. 

Because of differences between the experimental edge and the Gaussian pulse taken in part 1, 
and the limitations in filter size (it was not possible to obtain smaller filters), no conclusion can 
be drawn from the results, though a general agreement may be found. 
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Fig. 6 — Plots of the intensity distribution in the aerial image of an :T ’ 

unfiltered and two filtered edges (a = 1.5 and 2.5) 
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Fig. 7 — Normalized (to total area under unfiltered curve) intensity 
distribution of the image of two experimentally filtered edges 
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4. EDGE ENHANCEMENT THROUGH EXPOSURE ADDITION 


The next area to be investigated was that of actually applying filtering techniques in order to 
enhance an edge on a photographic emulsion. This edge may be enhanced in two different ways: 

(1) increase in contrast, and/or (2) increase in slope. 

Two methods for producing such enhancement using the image enhancement viewer have 
been considered. The first is called exposure addition. In this method, an unfiltered edge is ex- 
posed on an emulsion, after which an occluding filter is placed in the system and, without changing 
the position of either the edge or the film, a second exposure is made on the same frame. The 
second method is transmission multiplication. Here, two exposures are made on two separate 
frames; one exposure is unfiltered, the other is filtered to some degree. These two frames are 
then combined, using some sort of enlarging technique, in various proportions. Obviously, the 
difficulty in this method is the alignment of the two negatives, which is quite critical. This must 
be done with an accuracy that entails a great deal of technical ability. 

It was possible to do some work in the area of exposure addition, using the test edges from 
the previous area, but because of time limitation, it was impossible to do any transmission 
multiplication. 

An edge was exposed, using several occluding filters with varying exposure times, mixing 
the filtered and unfiltered parts. These edges were then microdensitometered on our Intectron 
instrument. For these traces, an effective 3 -micron slit was used. The results of these traces 
have been recorded in Table 1. 

Tables 1 and 2 are based upon exposures 1 to 19. The slopes (Table 1) and edge contrasts 
(Table 2) are those read, averaged, and weighted from several microdensitometer traces, varied 
for relative slit position with the edge. Included is the sensitometric curve relevant for the 
longer exposure times. Given also are the exposure times and the maximum densities (upper 
point through which the slope passes). As can be seen, some of these values are well up on the 
shoulder of the sensitometric curve (Fig. 8). All of the exposures were on 70-mm Panatomic-X 
film. A plot of the measured edge slopes against the ratio of filtered to unfiltered exposure time 
is plotted in Fig. 9. 

From this data it may be concluded that, with correct mixtures of filtered and unfiltered ex- 
posures, it is possible to equal or even slightly surpass the slope of the original edge. 
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Table 1 — Slope of Traced, Enhanced Edges 


Exposure 

Filter 

Average 

Unfiltered 

Filtered 

No. 

Diameter, in. 

Slope 

Time, sec 

Time, sec 

1 

— 

4.5 

1/2 

— 

2 

— 

3.5 

1 

— 

3 

0.015 

3.2 

1/2 

1 

4 

0.015 

2.5 

1/2 

2 

5 

0.015 

3.6 

1/2 

4 

6 

0.015 

4.3 

1 

1 

7 

0.015 

2.8 

1 

2 

8 

0.015 

2.4 

1 

4 

9 

0.025 

2.8 

1/2 

2 

10 

0.025 

4.0 

1/2 

4 

11 

0.025 

4.1 

1/2 

8 

12 

0.025 

2.8 

1 

4 

13 

0.025 

2.9 

1 

8 

14 

0.041 

2.8 

1/2 

4 

15 

0.041 

3.1 

1/2 

8 

16 

0.041 

5.0 

1/2 

16 

17 

0.041 

3.2 

1 

4 

18 

0.041 

2.0 

1 

8 

19 

0.041 

1.9 

1 
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Table 2 — 

Contrast of 

Traced, 



Enhanced Edges 


Exposure 

Minimum 

Maximum 


No. 

Density 

Density 

Contrast 

1 

0.60 

1.50 

0.90 

2 

0.95 

1.20 

0.85 

3 

0.85 

1.45 

0.60 

4 

1.35 

1.80 

0.45 

5 

1.55 

1.90 

0.35 

6 

1.15 

1.85 

0.70 

7 

1.50 

2.00 

0.50 

8 

1.40 

1.90 

0.50 

9 

0.80 

1.55 

0.75 

10 

1.20 

1.70 

0.50 

11 

1.00 

1.60 

0.60 

12 

1.10 

1.80 

0.70 

13 

1.35 

1.85 

0.50 

14 

0.85 

1.55 

0.70 

15 

0.75 

1.50 

0.75 

16 

0.90 

1.55 

0.65 

17 

0.95 

1.70 

0.75 

18 

1.30 

1.80 

0.50 

19 

1.65 

2.25 

0.60 
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Fig. 9 — Change of slope of an edge with a varying ratio of unfiltered time 
to filtered time using the technique of exposure addition 
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5. 


APPLICATION OF SPATIAL FILTERING TO THE 
OF AERIAL PHOTOGRAPHS 


EXAMINATION 


As part of the studies with the image enhancement viewer, it was considered desirable to 
examine some aerial negatives with various degrees of spatial filtering. The reason for this was 
that it is sometimes thought that details below the limit of detection can be brought above this 
threshold by suitable filtering. Also, it is conceivable that recognition of poorly defined objects 
might be improved by a controlled filtering which causes ringing at the edges. 

When the detail to be detected has a well defined spectrum, or a known pattern, spatial filter- 
ing has definite advantages. For example, a sinusoidal frequency of low modulation, recorded 
photographically, may be invisible against a grainy background, but by use of a spatial filter per- 
mitting only that frequency to pass, the grain noise is eliminated from the reconstituted image 
and the sinusoidal pattern may be seen. To be more precise, only those components of the grain 
frequency spectrum which lie in the passband of the filter remain to obscure the wanted image, 
and if the grain spectrum is flat, this gives a great improvement in signal -to-noise ratio. How- 
ever, this is an ideal case for the use of filtering, since the signal spectrum is concentrated at 
one point in the diffraction plane, while the noise spectrum is widely distributed across the plane. 
When the object is more complex— for example, a square wave instead of a sine wave — its spec- 
trum theoretically extends to infinity, and noise cannot be reduced without distorting the signal 
to some degree. If the significant feature of the object, however, is its repetitive nature and not 
the shape of the waves, the filtering can be as effective as for the sinusoidal pattern. For aerial 
scenes in general, little or no improvement can be expected from spatial filtering, because the 
details, consisting of sharp edged lines or similar, simple, nonrepetitive shapes, have very com- 
plex spectra. Filtering to reduce noise can only be effective if the objects of interest have sim- 
ple, known spectra, for which individual filters can be made in advance. 

To test these ideas under a realistic approximation to operational conditions, the following 
experiment was carried out with a setup which was being used in another program. Simulated 
aerial photographs were made by photographing a transparency made under controlled tone re- 
production conditions from a large scale aerial negative. The original negative scale was 1:2,000 
and the reduction factor, using a 6 -inch lens, was 20 to 1. Thus the copy negatives were at 
1:40,000 scale. Uniform fogging light from a beam splitter was introduced during the copying, to 
reduce the luminance scale of the transparency to approximately 5 to 1, simulating an aerial 
scene viewed from a high altitude. Negatives were made on SO-243 emulsion, varying the lens- 
film transfer function by changing the lens aperture. These lens -film transfer functions are 
shown in Fig. 10; for these values of spatial frequency cutoff, granularity is not obtrusive when 
looking for detail in the negatives. Thus a series of negatives were available in which the quality 
of definition showed a progressive deterioration of known amount. The intention was to apply 
various kinds of filtering, and if any improvement could be seen, to equate this to the quality of 
the next best negative. Thus any effect of the filtering could be approximately evaluated in terms 
of the definition of conventionally viewed photographs. 
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The original 1:40,000 scale negatives, being only 1 centimeter in diameter, were enlarged 
6 X , the resulting positive transparencies being more convenient for viewing in the spatial filter- 
ing apparatus. (It was known, from work in the other programs referred to above, that the en- 
larging process introduced no significant loss of quality even for the best negative.) 

After repeated trials of filtering, both high pass and low pass, it was concluded that no choice 
of frequency band enabled anything to be seen which could not be seen as well or better in normal 
viewing by incoherent light. Indeed, due to the disturbing effect of the ringing around small part- 
icles of dirt which had not been completely removed from the negatives or filtered from the fluid, 
the incoherent viewing was always preferred. In a few instances it was felt that the edges of 
isolated objects such as trucks might have been more clearly defined when using a high pass filter, 
but this effect was so small (if it existed) that a followup with more precise experimentation was 
not justified. 

This negative conclusion is in accordance with the theoretical expectations, as discussed 
earlier in this section. 

It should perhaps be stressed that these findings do not imply that coherent light techniques 
have no value for the detection of previously known patterns. 



Spatial Frequency, cycles per millimeter 
Fig. 10 — Four lens -film transfer functions using SO-243 film 
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A STUDY OF SPATIAL FILTERING USING 
THE IMAGE ENHANCEMENT VIEWER 


Introduction 

The work under the present phase of the contract consisted of two main 
parts. The first was the modification of the Image Enhancement Viewer and the 
testing and calibration of this instrument. The second was the continuance of 

research, in certain areas, which had been carried out by , while 

working under previous contracts for the Viewer. 

A description of the modified Viewer, along with the methods of testing 
and calibration is contained in the Operating Manual. Suffice it to say here 
that the Viewer was put in working order to the specifications listed in the 
manual. 

The following areas of research were investigated: 

A non-sharp, Gaussian pluse was theoretically passed through an optical 
coherent system; its low frequencies were filtered with a high-pass occluding 
filter, and the final images were plotted. 

A filmed non-sharp edge was experimentally filtered using the Viewer, 
and its image was traced with a scanning photometer. The results were curve 
fitted, normalized and plotted. 

Techniques (in this study, mainly exposure addition) for increasing 
edge contrast and sharpness using the Viewer were experimentally examined. 

The application of various filtering techniques using actual aerial 
photographs was studied. 


Approved For Release 2002/09/04 : CIA-RDP67B00820R0003001 10031-7 




Approved For Release 2002/09/045CIA-RDP67B00820R0003001 10031-7 


Finally, the use of occluding spatial filtering as a tool in the area 
of photographic mensuration was investigated. 

From the study of these five areas, a comprehensive view of the 
capabilities and limitations of the Image Enhancement Viewer and occluding 
spatial filtering may be obtained. 
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Theoretical Study of the Effect of Occluding Filters on an 
Unsharp (Gaussian) guise 

A full description of the Image Enhancement Viewer (Fig. 1) and its basic 
properties may be found in both the Operation Manual and the Itek Technical 
Report OD-61-6 . In operation it employs coherent illumination permitting a 
treatment that is linear in amplitude. 

Since all object films used in the system are immersed in an index- 
matching fluid (See Manual) there should be no shift in phase, and the object 
maybe considered solely in terms of amplitude. 

In order to obtain a general idea of the effect of occluding filtering 
on amplitude objects, a generalized case (considered Gaussian for mathematical 
simplicity) was considered theoretically. This pulse is transformed into the 
frequency plane and there filtered by a circular , centered occluding filter of 
generalized dimension. The resulting distribution is then again transformed 
into a filtered image, which has been computed and plotted for the range at 
representative values. 

Considering a non-sharp pulse of Gaussian form (Fig. 2) 


f x (x) 





2 

-cx 

e 


where f^(x) is the intensity distribution of a dark centered pulse in the 
object plane. Its amplitude spectrum in the frequency plane is given by 


CX) 



1/2 


- 1C0X , 

e dx 



1/2 


may 


then be rewritten in the following form 
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1/2 . 1/2 
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Expanding, we get 
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Letting the binomial coefficients be 


A(l) = 1 
2 

A (2) = 1 


2.4 

A(3) = 1.3 


A(4) = 1.3.5 

2 . 4 . 6. 8 

A(5) = 1.3. 5.7 

2.4.6.8.10 


2.4.6 



2 1/2 


n 

Then [l - (l - fl \ 1 - 1 

-£ A<n) 

A-M e — 

L ' 2 / -1 

n-1 

V T 2/ 

Putting this in the integral 
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Using the cosine transform, because of symmetry, and the aperture limits 


n cx> 


F(«) = 2 t 2 J 
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e cos o)X 


Integrating 
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This pulse may also have a light center in the form 


f 2 (x) = T 1 + [T 2 - Tj e 


(h - h) 


-cx 


The expansion and analysis using this pulse is very similar to that 
of f ^ (x) . 

We then obtain the image of f^(x), after passing it through a 
symmetric sharp cutoff occluding filter (Fig. 3) (centered in the frequency 
plane) whose characteristic is 


G(o>) = 1 




1 < “ <ex> 


= 0 


0 < 


The distribution in the image plane is given by the following transform 
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Again, for reasons of symmetry, we may use the cosine 


transform 
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Applying the Trigonometric identity 

sin(tox,) cos(toz) = 7 T sin(x. + z) w + sin(x_ - z) w 


and integrating, we get 
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I A(n) 

■y rt n=l 

Finally, the pulse image is given by 




Rewriting the image term in the following from where 
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Since c, the pulse variance is large, than x f W is large and the 
value of Sijx^V 2c + mj a approaches n/2. Under this condition, the 
image reduces to 


i^(a,m) = -V^ X! 
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1 " T, 
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At this point, it is necessary to compute the values of the integral 


cos(amt) dt 

Upon receiving these values, it is then possible to normalize the total 
energy of the filtered cases to the unfiltered case. 

The cases computed were those where "a" had the values of 1/4, 1/2, 

1 and 2 and n varied from 1 to 5 . m was taken in steps of .05 from 0 
to 3.0. was taken in steps of .01 from 0 to .99. T 2 was set at 1.0. 

The normalization procedure is the following 

h. 2 X*uml . I(a> m) 

2 

which is normalized to the total energy of the unfiltered pulse 

A = 1(0, m) where a = 0. Therefore, the normalization factor is 
P 

A 

N = T , P -r 
I(a ,m) 

The area of the filtered image, I(a,m) is approximated by Simpson's 
rule for normalization purposes. 

Representative results of the computations were plotted (Fig. 4 

a through h) . 

The parameters relating "k", the linear frequency, and "c", the 
pulse variance for "a's" varying from 0 to 2 are plotted in Fig. 5. 

This relationship is defined by 
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c 


2 2 
2n k 


2 

a 


For the entire range of edge contrast (Variation of T^), the 
smaller filter sizes give the expected results: increased normalized 
contrast for increased filter size. The larger filter sizes give rather 
unexpected and varying results for changing contrast. A full 
explanation for this phenomenon cannot be made without further study. 
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Experimental Filtering of a Non-Sharp Edge 

In order to substantiate the theoretical work, an experiment was 
designed with conditions approximating the analysis. 

Instead of a pulse, a filmed non-sharp edge was used as an object, being 
placed in the object plane of the Image Enhancement Viewer. Though this pulse 
had first been assumed to be of Gaussian shape, subsequent traces have revealed 
a number of "unGaussian" characteristics which could well affect the resulting 
filtered traces. 

Because an edge was used, and not a pulse, the filters placed in the 
object plane did not have to be circular (the frequency distribution being one 
dimensional. The occluding filters used were pieces of drill rod whose width 
was accurately known. These filters were then carefully centered in the 
frequency (or filter) plane. 

Using a scanning photometer whose slit width was minimized to the point 
where, at maximum photometer sensitivity, using the largest filter, 100/ o 
deflection on the recorder (5 mv range) could be reached. The edge was then 
traced ’with the photometer placed in the image plane with its slit parallel 
with the edge. First, the image of the edge with no occluding filter and the 
images of the edge with filters corresponding to a - 1.5 and a = 2.5 were traced 

The actual traces for the three cases are found in figure 6. (The two 
filtered cases are expanded by a factor of ten in the intensity scale: 

The resulting filtered cases were normalized to the area under the 
unfiltered curve in the same manner as used in the theoretical traces. The 
resulting curves are plotted in figure 7. Because the normalized cases are 
meant to be compared with the theoretical Gaussian, only the part of the edge 
corresponding to a pulse (center distance from the inflection point) is 

taken. 
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Because of differences between the experimental edge and the Gaussian 
pulse taken in Part I, and the limitations in filter size (it was not possible 
to obtain smaller filters), no conclusion may be drawn from the results, 
through a general agreement may be found. 


Approved For Release 2002/09/04 : CIA-RDP67B00820R0003001 10031-7 



Approved For Release 2002/09/04 : CIA-RDP67B00820R0003001 10031-7 


Edge Enhancement Through Exposure Addition 

The next area to be investigated was that of actually applying filtering 
techniques in order to "enhance" an edge on a photographic emulsion. This 
edge may be enhanced in two different ways; increase in contrast and/or 
increase in slope. 

At present, two methods for producing such enhancement using the Image 
Enhancement Viewer have been considered. The first is called exposure addition. 
In this method, an unfiltered edge is exposed on an emulsion, after which an 
occluding filter is placed in the system and without changing the position of 
either the edge or the film, a second exposure is made on the same frame. The 
second method is transmission multiplication. Here, two exposures are made on 
two separate frames; one exposure is unfiltered, the other, filtered to some 
degree. These two frames are then combined, using some sort of enlarging 
technique, in various proportions. Obviously, the difficulty in this method 
is the alignment of the two negatives which is quite critical. This must be 
done with an accuracy that would entail a great deal of technical ability. 

It was possible to do some work in the area of exposure addition, using 
the test edges from the previous area, but because of time limitation, it was 
impossible to do any transmission multiplication. 

The edge was exposed, using several occluding filters with varying 
exposure times, mixing the filtered and unfiltered parts. These edges were 
then microdensitometered on our Intectron instrument. For these traces, an 
effective 3 micron slit was used. The results of these traces have been 
recorded in Table #1. 
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Tables #1 and #2 refers to exposures #1-19. The slopes (Table #1) and 
edge contrasts (Table #2) are those read, averaged and weighted from several 
microdensitometer traces, varied for relative slit position with the edge. 
Included is the sensitometric curve relavent for the longer exposure times. 

Given also are the exposure times and the maximum densities (upper point through 
which the slope passes). As can be seen, some of these values are well up on 
the shoulder of the sensitometric curve (Table #3). All of the exposures were 
on Pan-X 70mm film. A plot of the measured edge slopes against the ratio of 
filtered to unfiltered exposure time is plotted in figure 8. 

From this data, it may be concluded that, with correct mixtures of 
filtered and unfiltered exposures, it is possible to equal, or even surpass, 
slightly, the slope of the original edge. 
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Application of Spatial Filtering to the Exami nation of Aerial Photographs 
As part of the studies with the Image Enhancement Viewer, it was 
considered desirable to examine some aerial negatives with various degrees of 
spatial filtering. The reason for this was that it is sometimes thought that 
details below the limit of detection can be brought above this threshold by 
suitable filtering. Also, it is conceivable that recognition of poorly defined 
objects might be improved by a controlled filtering which caused ringing at 

the edges. 

When the detail to be detected has a well-defined spectrum, or a known 
pattern, spatial filtering has definite advantages. For example, a sinusoidal 
frequency of low modulation, recorded photographically, may be invisible against 
a grainy background, but by use of a spatial filter permitting only that 
frequency to pass, the grain "noise" is eliminated from the reconstituted 
image and the sinusoidal pattern may be seen. To be more precise, only those 
components of the grain frequency spectrum which lie in the passband of the 
filter remain to obscure the wanted image, and if the grain spectrum is flat 
this gives a great improvement in signal-to-noise ratio. However, this is an 
ideal case for the use of filtering, since the signal spectrum is concentrated 
at one point in the diffraction plane while the noise spectrum is widely 
distributed across the plane. When the object is more complex, for example, a 
square wave instead of a sine wave, its spectrum theoretically extends to infinity, 
and noise cannot be reduced without distorting the signal in some degree. If the 
significant feature of the object, however, is its repetitive nature and not 
the shape of the waves, the filtering can be as effective as for the sinusoidal 
pattern. For aerial scenes in general, little or no improvement can be expected 
from spatial filtering, because the details, consisting of sharp edged lines 
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or similar simple non-repet itive shapes, have very complex spectra. Filtering 
to reduce noise can only be effective if the objects of interest have simple known 
spectra, for which individual filters can be made in advance. 

To test these ideas under a realistic approximation to practical conditions, 
the following experimentation was carried out, with a set-up which was being 
used in another program. Simulated aerial photographs were made by photographing 
a transparency made under controlled tone-reproduction conditions from a large- 
scale aerial negative. The original negative scale was 1:2000 and the reduction 
factor, using a 6 inch lens, was 20 to 1, thus the copy negative were at 1:40,000 
scale. Uniform fogging light from a beam-splitter was introduced during the 
copying, to reduce the luminance scale of the transparency to approximately 5 to 
1, simulating an aerial scene viewed from a high altitude. Negatives were made 
on SO 243 emulsion, varying the lens/film transfer function by changing the lens 
aperture. These lens/film transfer functions are shown in figure 9; for these 
values of spatial frequency cut-off, granularity is not obtrusive when looking 
for detail in the negatives. Thus a series of negatives were available in which 
the quality of definition showed a progressive deterioration of known amount. The 
intention was to apply various kinds of filtering, and if any improvement could 
be seen, to equate this to the quality of the next best negative. Thus any effect 
of the filtering could be approximately evaluated in terms of the definition of 
conventionally-viewed photographs. 

The original 1:40,000 scale negatives, being only 1 cm in diameter, were 
enlarged X6, the resulting positive transparencies being more convenient for 
viewing in the spatial-filtering apparatus. (It was known, from work in the 
other programs referred to above, that the enlarging process introduced no 
significant loss of quality even for the best negative.) 
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After repeated trials of filtering, both high-pass & low-pass, it was 
concluded that no choice of frequency band enabled anything to be seen which 
could not be seen as well or better in normal viewing by incoherent light. Indeed, 
due to the disturbing effect of the ringing around small particles of dirt which 
had not been completely removed from the negatives or filtered from the fluid, 
the incoherent viewing was always preferred. In a few instances it was felt 
that the edges of isolated objects such as trucks might have been more clearly 
defined when using a high-pass filter, but this effect was so small, (if it 
existed) that a follow-up with more precise experimentation was not justified. 

This negative conclusion is in accordance with the theoretical expectations, 

as discussed earlier in this section. 

It should perhaps be stressed that these findings do not imply that 
coherent light techniques have no value for the detection of previously-known 

patterns . 
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CHARACTERISTICS OF TRACED ENHANCED EDGES 
Table #1 


posur e 

Filter 

Average 

Unfiltered 

Filtered 

No . 

(IN) 

Slope 

Time 

Time 

i 


4.5 

1/2 

- 

2 

- 

3.5 

1 

- 

3 

.015 

3.2 

1/2 

1 

4 

.015 

2.5 

1/2 

2 

5 

.015 

3.6 

1/2 

4 

6 

.015 

4.3 

1 

1 

7 

.015 

2.8 

1 

2 

8 

.015 

2.4 

1 

4 

9 

.025 

2.8 

1/2 

2 

10 

.025 

4.0 

1/2 

4 

11 

.025 

4.1 

1/2 

8 

12 

.025 

2.8 

1 

4 

13 

.025 

2.9 

1 

8 

14 

.041 

2.8 

1/2 

4 

15 

.041 

3.1 

1/2 

8 

16 

.041 

5.0 

1/2 

16 

17 

.041 

3.2 

1 

4 

18 

.041 

2.0 

1 

8 

19 

.041 

1.9 

1 

16 
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Table #2 


Exposure # 


Min. 


Max. 


Contrast 


1 

.60 

2 

.95 

3 

.85 

4 

1.35 

5 

1.55 

6 

1.15 

7 

1.50 

8 

1.40 

9 

.80 

10 

1.20 

11 

1.00 

12 

1.10 

13 

1.35 

14 

.85 

15 

.75 

16 

.90 

17 

.95 

18 

1.30 

19 

1.65 


1.50 

.90 

1.20 

.85 

1.45 

.60 

1.80 

.45 

1.90 

.35 

1.85 

.70 

2.00 

.50 

1.90 

.50 

1.55 

.75 

1.70 

.50 

1.60 

.60 

1.80 

.70 

1.85 

.50 

1.55 

.70 

1.50 

.75 

1.55 

.65 

1.70 

.75 

1.80 

.50 

2.25 

.60 
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Major Features and Dimensions of 
Image Enhancement Viewer 

(Approximate overall system length « 8 feet) 
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GAUSSIAN SPATIAL PULSE WITH OCCLUDING SPATIAL FILTER 
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Fig 4a 


NORMALIZED IMAGE INTENSITY 
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NORMALIZED IMAGE INTENSITY 



FORMALIZED IMAGE INTENSITY 




NORMALIZED IMAGE INTENSU Y 
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NORMALIZED IMAGE INDENSITY 



uNFILTERED TIME/FILTERED TIME 
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Modulation Transfer 
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Spatial Frequency, cycles per millimeter 


Four lens -film transfer functions: SO-243 


Fig. 9 
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